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The first-principles density functional molecular dynamics simulations have been carried out to 
investigate the geometric, the electronic, and the finite temperature properties of pure Li clusters 
(Liio, Lii2) and Al-doped Li clusters (LiioAl, LiioAfj). We find that addition of two Al impurities 
in Liio results in a substantial structural change, while the addition of one Al impurity causes a 
rearrangement of atoms. Introduction of Al-impurities in Liio establishes a polar bond between 
Li and nearby Al atom(s), leading to a multicentered bonding, which weakens the Li-Li metallic 
bonds in the system. These weakened Li-Li bonds lead to a premelting feature to occur at lower 
temperatures in Al-doped clusters. In LiioAi2, Al atoms also form a weak covalent bond, resulting 
into their dimer like behavior. This causes Al atoms not to 'melt' till 800 K, in contrast to the Li 
atoms which show a complete diffusive behavior above 400 K. Thus, although one Al impurity in 
Liio cluster does not change its melting characteristics significantly, two impurities results in 'surface 
melting' of Li atoms whose motions are confined around Al dimer. 

PACS numbers: 61.46.Bc, 36.40.-c, 36.40.Cg, 36.40.Ei 



I. INTRODUCTION 

There have been considerable experimental and the- 
oretical studies to understand the physical and chem- 
ical properties of cluster which include structural and 
electronic properties, the nature of bonding, thermody- 
namics^ and spectroscopic properties. Since most of 
these studies have been carried out on the homogeneous 
clusters^ the physics of mixed clusters remains less ex- 
plored^ It is well known that the dilute impurities alter 
the electronic structure and geometries of the bulk sys- 
tem. A similar phenomenon is also observed in impurity- 
doped clusters 4£ 

Cheng et al. have investigated the energetics and the 
electronic structure of small Li-Al clusters £ They sug- 
gested a special role of AlLis unit as a building block for 



clusters of assembled materials, e.g. Al„Li5„ (n > 1). 
This idea was further explored by Akola and Manni- 
neni They investigated small Li-rich AljvLisjv (N=l- 
6,10) clusters using first-principle calculations. They re- 
ported that Al ions form a compact inner core embedded 
in Li atoms. However, they did not find AlLis to be a 
favorable candidate as a building block for larger clus- 
ters. Further, they observed a significant charge transfer 
from Li to nearby Al atoms, strengthening ionic bonds 
between Li and Al, as well as a formation of Al-Al cova- 
lent bonds. Similar findings on LiioAbj cluster have been 
reported by Kumar. 8 

Another issue of considerable interest is the finite tem- 
perature behavior of homogeneous as well as impurity- 
doped clusters. Joshi et al. have investigated the finite 
temperature behavior of impurity-doped cluster, LigSn. 5 
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Their work indicates that the addition of one impurity re- 
sults in lowering of melting temperature by about 125 K. 
A similar observation has made be Aguado et al. for 
the case of LiNas4 and CsNa54 clusters^ Recently, an in- 
teresting study on the effect of a single impurity in the 
icosahedral clusters of silver has been reported by Mot- 
tct et alw& In contrast to the previous studies;^ they 
showed that a single impurity of Ni or Cu can lead to an 
increase in the melting temperature of the host. A recent 
study by Zorriassatein et al. on the melting of Si^Ti also 
reveals that a single impurity like Ti can change the finite 
temperature behavior of the host cluster significantly. 11 
Therefore, it is of considerable interest to investigate the 
impurity induced effects in Li clusters. In the present 
work, we have carried out the first-principles density 
functional molecular dynamics simulations on LiioAl and 
LiioAl2 clusters. The results have been compared with 
those of pure Li clusters. In particular, we have investi- 
gated the equilibrium geometries, the nature of bonding, 
and the finite temperature behavior of these clusters. In 
sec. II, we briefly describe the computational procedure. 
We present our results and discussion in sec. III. A brief 
summary of results is given in sec. IV. 

II. COMPUTATIONAL DETAILS 

We have employed Bohn-Oppenheimer molecular dy- 
namics 12 using Vanderbilt's ultrasoft pseudopotentialsJ^ 
within the local density approximation, as implemented 
in the VASP package. 14 - A cubic supercell of length 20 A 
with energy cutoff of 9.5 Ry was used for the total-energy 
convergence. We have obtained the lowest energy struc- 
ture and other equilibrium geometries by a two step pro- 
cess. In the first step, starting from random configu- 
rations, the clusters are heated to a few representative 
temperatures, below and above expected melting points. 
The clusters are maintained at these temperatures for 



at least 60 ps. Then, resulting trajectories are used to 
choose several initial configurations for geometry opti- 
mizations. In this way, we have obtained various equi- 
librium geometries. The geometries are considered to 
be converged when the force on each ion is less than 
0.005 eV/A with a convergence in the total energy to 
be the order of 10 eV. To investigate the nature of the 
bonding, we have examined the total charge density and 
the molecular orbitals (MO's). 

To examine the finite temperature behavior of clusters, 
molecular dynamics simulations have been carried out at 
16 temperatures for Liio, and at 13 temperatures for Lii2 
and LiioAl m (m=l,2) within the range of 100 K < T < 
800 K. The simulation time for each temperature is at 
least 150 ps. We have discarded the first 30 ps for each 
temperature to allow the system to be thermalized. The 
resulting trajectories have been used to calculate stan- 
dard thermodynamic indicators as well as the ionic spe- 
cific heat via multiple histogram technique. The details 
can be found in Ref. 

III. RESULTS AND DISCUSSION 
A. The Geometries and the Electronic Structures 

The lowest energy structures and various isomers of 
Liio, Lii2, LiioAl, and LiioA^ are shown in Figs.^and 
12 The ground state (GS) geometry of Liio consists of 
two interconnected pentagonal rings with planes perpen- 
dicular to each other (Fig. d-a(i)), which agrees with 
that found by Fournier et alw& as well as with our earlier 
study^ The structure with an atom at the center (Fig. 
HJ-a(ii)) is nearly degenerate with the ground state (Fig. 
IU-a(i)). This is the structure reported by Jones et al. as 
their ground stated A high energy structure shows an 
interconnection of dodecahedron and octahedron (Fig.EJ- 
a(iii)). We have found five nearly degenerate equilibrium 
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FIG. 1: The ground-state geometry and some isomers of Liio 
and Lii2. The label (i) represents the ground-state geometry. 
The energy difference AE is given in eV with respect to the 
ground-state energy. 

geometries of Lii2 within the energy range of 0.007 eV. 
We examined the stability of these structures by vibra- 
tional analysis. Three of these are shown in Figs. ID-b(i) 
~n~b(iii), where Fig.2]-b(ii) have been reported as the 
lowest-energy structure by Fournier et alw& All the struc- 
tures consist of two units: a pentagonal bipyramid and a 
distorted octahedron connecting to each other with dif- 
ferent angles. As we shall see, the existence of these 
nearly degenerate structures have a bearing on the shape 
of the specific heat curve at low temperatures. In high 
energy structures, the octahedron is destroyed first (Figs. 
[Tl-b(iv) andn~b(v)) an d then the destruction of pentag- 
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FIG. 2: The ground-state geometry and some isomers of 
LiioAl and LiioAla. The label (i) represents the ground-state 
geometry. The blue circle represents the Li atoms and the 
yellow circle represents the Al atoms. The energy difference 
AE is given in eV with respect to the ground-state energy. 

onal bipyramid follows (Fig.[IJ-b(vi)). 

Now, we discuss the geometries of impurity-doped sys- 
tems. So far as Liio is concerned, a single Al atom re- 
places one of the Li atoms in one unit of the pentagonal 
bipyramids in the GS geometry of LiioAl (Fig. EJ-c(i)). 
This causes the Li atoms in another unit to rearrange in 
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the form of antiprism, with Al atom at the center. As 
shown in Figs. EJ-c(ii) to|2-c(v), the low energy geome- 
tries are dominated by the presence of pentagonal bipyra- 
mid, whereas the high energy structures are dominated 
by antiprism. The addition of one more Al atom in LiioAl 
changes the lowest energy structure significantly. The 
equilibrium structures of LiioAL) do not show the pen- 
tagonal bipyramidal structure any more. Instead, there 
is an octahedron consisting of four Li atoms and two 
Al atoms with Li atoms forming a central plane. This 
core octahedron is common to all the isomers. Different 
isomers represent different ways of capping this core by 
the remaining Li atoms. Symmetric arrangements of Li 
atoms give rise to low energy isomers (Fig. G]-d(i) ~ 0- 
d(iii)). The lowest energy structure (Fig. El-d(i)) is the 
most symmetric structure among all the clusters studied. 
This is because two Al atoms share Li atoms equally to 
fill their unoccupied p-orbitals. The nearly degenerate 
structure shown in Fig. |2-d(ii) has been found as the 
lowest energy structure by Cheng et The stability 
of two nearly degenerate structures has been verified by 
carrying out a vibrational analysis. The structures with 
broken symmetry have higher energies (Figs. El-d(iv) ~ 
E-d(vi)). 

We used the deformation parameter Edef to examine 
the shape of the clusters. Edef is defined as, 

2Qi 

£def ~ Q2 + Q3 
where Q\ > Q2 > Q3 are eigenvalues of the quadrupole 
tensor 

Qij = ^ RaRij 
1 

with Rn being i th coordinate of ion / relative to the 
center of mass of the cluster. A spherical system has 
Edef = 1 (<3i = Q2 = Q3), while Edef > 1 indicates a 
deformation. The calculated Edef for the GS geometry of 
the clusters are shown in Fig. |3{a). It can be seen that 
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FIG. 3: The properties of the ground state geometry, (a) The 
deformation parameter (edef)', (b) The binding energy Ef, in 
eV/atom. 

doping reduces the deformation considerably. We also 
show the binding energies of these clusters in Fig. Efb). 
Clearly the increase in the binding energy by adding Al 
atom (0.2 eV/atom) is much higher than that by adding 
Li atom (0.04 eV/atom). Thus, we conclude that the 
increase of the binding energy in the mixed clusters is 
mainly due to the formation of strong Li-Al bonds. 

We examine the change in the nature of the bonding 
due to the presence of Al atoms via the total electron 
charge density, and the molecular orbitals (MO's). In 
addition, we have also calculated the difference charge 
density between electron charge densities of mixed clus- 
ter LiioAl m (to=1,2) and separated units of Liio and Al m 
(m=l,2), by keeping the atomic positions the same as 
those in the mixed cluster. We first examine the charge 
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FIG. 4: The isovalued surfaces of the total charge density 
at 1/3 of its maximum values, where the maximum charge 
density of LiioAl is 0.233 and that of LiioA^ is 0.239. 

distribution in AI2 dimcr and Li-Al diatomic cluster. Ev- 
idently, there is a covalent bond in AI2 dimer as expected, 
while Li-Al diatomic cluster shows a polar bond. This 
can be explained on the basis of electronegativity differ- 
ence between Li(0.98) and Al(1.61), which is not large 
enough to establish an ionic bond. We note that the 
bond lengths of these clusters are 2.68 A for Al-Al and 
2.88 A for Li-Al. The isovalued surface of total charge 
density in Li^Al and LiioA^ are shown in Figs. EJa) 
and^Jb), respectively. In Figs. EKa) toEJd), we show 
the isosurface of difference charge density for these Al- 
doped clusters. When we add one Al atom to the Liio 
cluster, the total charge distribution shows a mixed char- 
acter of localization near Al atom and derealization on 
the pentagonal bipyramid in Li^Al (Fig. Ufa)). It can 
be seen that the charge, localized around Al atom in the 
antiprism part, is nearly spherical. This is due to a sig- 
nificant charge transfer from nearby Li atoms to the Al 
atom to fill its unoccupied p-orbitals However, there is 
an insignificant charge transfer from the Li atoms to the 
Al atom, forming the pentagonal bipyramidal unit. This 
leads to a delocalized charge distribution in this unit. 
This observation is confirmed by the examination of Figs. 
EJa) andj^b), depicting the constant density contour of 
the difference charge density. It may be noted that the 
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FIG. 5: The difference charge density between the mixed clus- 
ter and separated units of Li and Al clusters at the same 
positions, (a) and (c) show the region where the charge is 
gained and (b) and (d) show the region where the charge is 
lost. The figures are shown the isosurfaces of difference charge 
density at the value of charge gain of 0.017 and loss of 0.008, 
respectively. 

values of the maxima in the difference charge density are 
0.042 and 0.051, and those of the minima are -0.037 and 
-0.084 in LiioAl and LiioA^, respectively. Figs.OJa) to 
I3d) show the isosurfaces for the value of 0.017 for the 
charge gain, and -0.008 for the charge loss. It can be seen 
that a ring shape region in between Li and Al (around 
Al) gains charge in LiioAl (Fig. Eta)). Further, most 
of the charge is lost by Li atoms nearer to the Al atom 
(Fig-Etb)). This charge transfer from Li to Al causes Li 
atoms to be positively charged. These positively charged 
Li atoms pull the charge distribution around Al atom 
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to polarize it. This leads to a multicentered bonding 
between Al atom and nearby Li atoms. The addition 
of one more Al atom changes the charge density distri- 
bution from a mixed localized and delocalizcd one to a 
mainly localized one. In LiioAl2, the charge distribution 
(Fig- Elk)) is mainly around two Al atoms polarized by 
nearby Li atoms. The difference charge density shows 
a very symmetric charge gain region in between Li and 
Al atoms (Fig. EJc)) and loss of charge at each atomic 
site (Fig.Efd)). These difference in the bonding between 
LiioAl and LiioAl2 are also seen in molecular orbitals 
(figures not shown). For example, in LiioAl three high- 
est occupied molecular orbital (HOMO) show a dclocal- 
ized charge distribution, while only the HOMO shows a 
derealization in LiioA^. As noted earlier, since elec- 
tronegativity difference of Li and Al is not very high, the 
Li atoms donate their charge to Al atom(s) partially. We 
expect this charge transfer to result in weakening of Li- 
Li metallic bond in the system. This feature has also 
been noted in LigSn cluster £ Further, this partial charge 
transfer also results in the sharing of charge between two 
Al atoms to fill their unoccupied p-orbitals in LiioA^. 
As we shall see, this difference in the bonding leads to a 
different thermodynamic behavior in two clusters. 

B. Thermodynamics 

1. Pure Li clusters 

To investigate the melting behavior of the clusters, we 
calculate the canonical specific heat using multiple his- 
togram technique. The calculated ionic specific heats for 
all the clusters are shown in Fig. El It is well known 
that a small cluster exibits a broad melting transition. 
Our observatios are consistent with this observation. In 
addition, in Liio the ionic specific heat also shows a re- 
markable premelting feature between 150 K and 225 K, 
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FIG. 6: The normalized specific-heat as a function of temper- 
ature. Co = (3iV — 9/2)fcs is the zero-temperature classical 
limit of the rotational plus vibrational canonical specific-heat. 



whose maximum value is close to that of the main peak 
around 575 K. The examination of the ionic motion re- 
veals that the two structures shown in Figs. H~a(i) and 
n~a(ii) are observed in the temperature range of 150 K 
~225 K. The system visits these two isomers in the time 
spent of a few pico seconds at 175 K, indicating that 
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this premelting feature is due to the isomerization. At 
higher temperature of 260 K, we observe the structure 
shown in Fig. |IJ-a(iii). At still higher temperatures, the 
cluster visits these three isomers frequently. It is diffi- 
cult to identify the melting temperature for this cluster 
because liquidlike behavior develops over a wide range of 
temperature (300 K to 700 K). The root-mean-square 
bond length fluctuation (S rms ) is another indicator for 
studying a melting transition. According to the Linde- 
mann criteria (for bulk) solid-liquid transition is signified 
when the value of S rms exceeds 0.1. However, it is gener- 
ally observed that for clusters a liquidlike behavior is seen 
when 5 rms exceeds 0.25^0.3. As shown in Fig. S rms 
increases in two steps. At 175 K it exceeds the value of 
0.1, which corresponds to maximum specific heat of the 
shoulder. However, as discussed earlier, the cluster is not 
in the liquidlike state. S rms increases from 200 K again 
until it saturates about 0.3 at 575 K where we observe 
the maximum value of the specific heat. 

The specific heat of Lii2 shows a weak shoulder around 
125 K before it increases gradually. This shoulder is re- 
lated to the existence of five degenerate isomers. At low 
temperatures up to 150 K, the motion is dominated by 
vibrations and the cluster visits the low energy isomers. 
As temperature increases, the liquidlike behavior evolves 
over a very wide range of temperature. At 175 K, we ob- 
serve high energy isomers where the octahedral structure 
is absent (two of them are shown in Figs.^b(iv) and^ 
b(v)). With further increase of temperature to 375 K, we 
observe that the pentagonal bipyramidal structure is de- 
stroyed (one of these structures is shown in Fig.^-b(vi)). 
It can be seen that 5 rms also increases continuously from 
125 K to 500 K (Fig. Q). It is interesting to note that 
even though the isomerization is seen in both the clusters 
Liio and Lii2 , evidently the isomerization in Liio leads to 
a prominent shoulder in its specific heat. This is related 



to the nature of the process by which their ground state 
visits the isomer(s). In the case of Liio, there is a barrier 
at about 175 K. The system overcomes this barrier with 
sudden increase in the accessible density of states. For 
Lii2, since the different isomers differ only in the angle 
between two constituent units, the ground state visits 
them almost continuously, aided by the vibrational mo- 
tion. This results in different specific heat curves. 

2. Thermodynamics of Al-doped Li clusters 

In can be seen from Fig.|Blthat LiioAl shows a premelt- 
ing feature in the low temperature region around 150 K. 
This may be caused by a weakening of the Li-Li metal- 
lic bond due to the charge transfer from Li atoms to the 
Al atom in antiprism unit. This effect has also been ob- 
served in LigSn£ The Li atoms in this unit break the 
Li-Li bonds at low temperature and one of them moves 
towards the pentagonal bipyramidal unit. This is the way 
the cluster visits the other low energy isomers shown in 
Figs. I2]-c(ii) and 0-c(iii). For the temperatures above 
175 K, we observe the motion of Li atoms around the Al 
atom without destroying the overall shape of the struc- 
ture. However, the four Li atoms consisting of pentagonal 
bipyramidal unit do not interchange their position with 
others till about 575 K. At this temperature, the sys- 
tem shows a diffusive liquidlike behavior. We calculate 
the root-mean-square bond length fluctuation (S rms ) for 
Li-Li, Li-Al, and Al-Al separately, to see the difference 
between their fluctuations. 5 rrns is defined as, 

where Rij is the distance between i th and j th ions with i 
and j for relevant bonds, and N is the number of bonds. 

denotes a time average over the entire trajectory. It 
can be seen in Fig. that in the low temperature region 
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FIG. 7: The root-mean-square bond length fluctuation 
(firms) of Li-Li, Li-Al, Al-Al as a function of temperature. 

5 rrns of Li-Al is much lower than that of Li-Li. This is 
because of relatively stronger Li-Al bond. 

Interestingly, the specific heat of Li 10 Al 2 at low tem- 
perature also shows a shoulder which is very similar to 
the specific heat of Li 10 . After this premelting feature, 
the specific heat remains flat. This is due to the presence 
of Al dimer in the cluster. When we compare its molecu- 



lar orbitals with those of pure Al dimer, they looked alike 
even though their bond lengths are different. The bond 
length of Al-Al in LiioAi2 is 2.97A, while that of pure 
Al dimer is 2. 68 A. Thus, there is a weak covalent bond 
between Al atoms. This covalent bond between them re- 
stricts the motion of Al atoms. We observe that Al-Al 
bond does not break at least up to the temperature of 
about 800 K. As shown in Fig. [7] 5 rms of Li-Li and Li-Al 
saturate to a value of about 2.5 at 500 K, while that of 
Al-Al reaches a value of 0.1 at 800 K. The saturation 
values of 5 rms for Li-Li and Li-Al are lower than corre- 
sponding S rms of other clusters, indicating the motion of 
atoms in LiioAi2 is more restricted than other systems 
at comparable temperature. 

IV. SUMMARY 

We have employed ab initio molecular dynamics to 
study the equilibrium geometries, the electronic struc- 
ture, and the finite temperature properties of pure 
lithium clusters Liio and Lii2 , and Al-doped lithium clus- 
ters LiioAl and LiioAk>. We find that there is a substan- 
tial structural change upon doping with two Al atoms, 
while the addition of one Al atom results in a rearrange- 
ment. The analysis of the total charge density and the 
molecular orbitals reveal that there is a partial charge 
transfer from Li atoms to Al atom(s) in Al-doped clus- 
ters, forming a polar bond between them. This leads to 
a multicentered bonding and weakens the Li-Li metallic 
bonds in these clusters. In LiioA^, Al atoms also form 
a weak covalent bond. These changes in the nature of 
bonding upon doping affect the finite temperature prop- 
erties of pure host cluster. We observe that the presence 
of dimer - like Al atoms with weak covalent bond confines 
the motion of Li atoms around them. Thus, a substi- 
tution of two Li atoms by Al atoms in Lii2 leads to a 
surface melting only, showing a continuous phase change 
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